\(\int \frac {1-x^4}{1+b x^4+x^8} \, dx\) [20]
Optimal result
Integrand size = 20, antiderivative size = 511 \[
\int \frac {1-x^4}{1+b x^4+x^8} \, dx=-\frac {\sqrt {2+b} \arctan \left (\frac {\sqrt {2-\sqrt {2-b}}-2 x}{\sqrt {2+\sqrt {2-b}}}\right )}{4 \sqrt {2-\sqrt {2-b}} \sqrt {2-b}}+\frac {\sqrt {2+b} \arctan \left (\frac {\sqrt {2+\sqrt {2-b}}-2 x}{\sqrt {2-\sqrt {2-b}}}\right )}{4 \sqrt {2+\sqrt {2-b}} \sqrt {2-b}}+\frac {\sqrt {2+b} \arctan \left (\frac {\sqrt {2-\sqrt {2-b}}+2 x}{\sqrt {2+\sqrt {2-b}}}\right )}{4 \sqrt {2-\sqrt {2-b}} \sqrt {2-b}}-\frac {\sqrt {2+b} \arctan \left (\frac {\sqrt {2+\sqrt {2-b}}+2 x}{\sqrt {2-\sqrt {2-b}}}\right )}{4 \sqrt {2+\sqrt {2-b}} \sqrt {2-b}}+\frac {\sqrt {2-\sqrt {2-b}} \log \left (1-\sqrt {2-\sqrt {2-b}} x+x^2\right )}{8 \sqrt {2-b}}-\frac {\sqrt {2-\sqrt {2-b}} \log \left (1+\sqrt {2-\sqrt {2-b}} x+x^2\right )}{8 \sqrt {2-b}}-\frac {\sqrt {2+\sqrt {2-b}} \log \left (1-\sqrt {2+\sqrt {2-b}} x+x^2\right )}{8 \sqrt {2-b}}+\frac {\sqrt {2+\sqrt {2-b}} \log \left (1+\sqrt {2+\sqrt {2-b}} x+x^2\right )}{8 \sqrt {2-b}}
\]
[Out]
-1/4*arctan((-2*x+(2-(2-b)^(1/2))^(1/2))/(2+(2-b)^(1/2))^(1/2))*(2+b)^(1/2)/(2-b)^(1/2)/(2-(2-b)^(1/2))^(1/2)+
1/4*arctan((2*x+(2-(2-b)^(1/2))^(1/2))/(2+(2-b)^(1/2))^(1/2))*(2+b)^(1/2)/(2-b)^(1/2)/(2-(2-b)^(1/2))^(1/2)+1/
8*ln(1+x^2-x*(2-(2-b)^(1/2))^(1/2))*(2-(2-b)^(1/2))^(1/2)/(2-b)^(1/2)-1/8*ln(1+x^2+x*(2-(2-b)^(1/2))^(1/2))*(2
-(2-b)^(1/2))^(1/2)/(2-b)^(1/2)+1/4*arctan((-2*x+(2+(2-b)^(1/2))^(1/2))/(2-(2-b)^(1/2))^(1/2))*(2+b)^(1/2)/(2-
b)^(1/2)/(2+(2-b)^(1/2))^(1/2)-1/4*arctan((2*x+(2+(2-b)^(1/2))^(1/2))/(2-(2-b)^(1/2))^(1/2))*(2+b)^(1/2)/(2-b)
^(1/2)/(2+(2-b)^(1/2))^(1/2)-1/8*ln(1+x^2-x*(2+(2-b)^(1/2))^(1/2))*(2+(2-b)^(1/2))^(1/2)/(2-b)^(1/2)+1/8*ln(1+
x^2+x*(2+(2-b)^(1/2))^(1/2))*(2+(2-b)^(1/2))^(1/2)/(2-b)^(1/2)
Rubi [A] (verified)
Time = 0.25 (sec) , antiderivative size = 511, normalized size of antiderivative = 1.00, number of
steps used = 19, number of rules used = 6, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.300, Rules used = {1435, 1183, 648, 632, 210,
642} \[
\int \frac {1-x^4}{1+b x^4+x^8} \, dx=-\frac {\sqrt {b+2} \arctan \left (\frac {\sqrt {2-\sqrt {2-b}}-2 x}{\sqrt {\sqrt {2-b}+2}}\right )}{4 \sqrt {2-\sqrt {2-b}} \sqrt {2-b}}+\frac {\sqrt {b+2} \arctan \left (\frac {\sqrt {\sqrt {2-b}+2}-2 x}{\sqrt {2-\sqrt {2-b}}}\right )}{4 \sqrt {\sqrt {2-b}+2} \sqrt {2-b}}+\frac {\sqrt {b+2} \arctan \left (\frac {\sqrt {2-\sqrt {2-b}}+2 x}{\sqrt {\sqrt {2-b}+2}}\right )}{4 \sqrt {2-\sqrt {2-b}} \sqrt {2-b}}-\frac {\sqrt {b+2} \arctan \left (\frac {\sqrt {\sqrt {2-b}+2}+2 x}{\sqrt {2-\sqrt {2-b}}}\right )}{4 \sqrt {\sqrt {2-b}+2} \sqrt {2-b}}+\frac {\sqrt {2-\sqrt {2-b}} \log \left (-\sqrt {2-\sqrt {2-b}} x+x^2+1\right )}{8 \sqrt {2-b}}-\frac {\sqrt {2-\sqrt {2-b}} \log \left (\sqrt {2-\sqrt {2-b}} x+x^2+1\right )}{8 \sqrt {2-b}}-\frac {\sqrt {\sqrt {2-b}+2} \log \left (-\sqrt {\sqrt {2-b}+2} x+x^2+1\right )}{8 \sqrt {2-b}}+\frac {\sqrt {\sqrt {2-b}+2} \log \left (\sqrt {\sqrt {2-b}+2} x+x^2+1\right )}{8 \sqrt {2-b}}
\]
[In]
Int[(1 - x^4)/(1 + b*x^4 + x^8),x]
[Out]
-1/4*(Sqrt[2 + b]*ArcTan[(Sqrt[2 - Sqrt[2 - b]] - 2*x)/Sqrt[2 + Sqrt[2 - b]]])/(Sqrt[2 - Sqrt[2 - b]]*Sqrt[2 -
b]) + (Sqrt[2 + b]*ArcTan[(Sqrt[2 + Sqrt[2 - b]] - 2*x)/Sqrt[2 - Sqrt[2 - b]]])/(4*Sqrt[2 + Sqrt[2 - b]]*Sqrt
[2 - b]) + (Sqrt[2 + b]*ArcTan[(Sqrt[2 - Sqrt[2 - b]] + 2*x)/Sqrt[2 + Sqrt[2 - b]]])/(4*Sqrt[2 - Sqrt[2 - b]]*
Sqrt[2 - b]) - (Sqrt[2 + b]*ArcTan[(Sqrt[2 + Sqrt[2 - b]] + 2*x)/Sqrt[2 - Sqrt[2 - b]]])/(4*Sqrt[2 + Sqrt[2 -
b]]*Sqrt[2 - b]) + (Sqrt[2 - Sqrt[2 - b]]*Log[1 - Sqrt[2 - Sqrt[2 - b]]*x + x^2])/(8*Sqrt[2 - b]) - (Sqrt[2 -
Sqrt[2 - b]]*Log[1 + Sqrt[2 - Sqrt[2 - b]]*x + x^2])/(8*Sqrt[2 - b]) - (Sqrt[2 + Sqrt[2 - b]]*Log[1 - Sqrt[2 +
Sqrt[2 - b]]*x + x^2])/(8*Sqrt[2 - b]) + (Sqrt[2 + Sqrt[2 - b]]*Log[1 + Sqrt[2 + Sqrt[2 - b]]*x + x^2])/(8*Sq
rt[2 - b])
Rule 210
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(-(Rt[-a, 2]*Rt[-b, 2])^(-1))*ArcTan[Rt[-b, 2]*(x/Rt[-a, 2])
], x] /; FreeQ[{a, b}, x] && PosQ[a/b] && (LtQ[a, 0] || LtQ[b, 0])
Rule 632
Int[((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(-1), x_Symbol] :> Dist[-2, Subst[Int[1/Simp[b^2 - 4*a*c - x^2, x], x]
, x, b + 2*c*x], x] /; FreeQ[{a, b, c}, x] && NeQ[b^2 - 4*a*c, 0]
Rule 642
Int[((d_) + (e_.)*(x_))/((a_.) + (b_.)*(x_) + (c_.)*(x_)^2), x_Symbol] :> Simp[d*(Log[RemoveContent[a + b*x +
c*x^2, x]]/b), x] /; FreeQ[{a, b, c, d, e}, x] && EqQ[2*c*d - b*e, 0]
Rule 648
Int[((d_.) + (e_.)*(x_))/((a_) + (b_.)*(x_) + (c_.)*(x_)^2), x_Symbol] :> Dist[(2*c*d - b*e)/(2*c), Int[1/(a +
b*x + c*x^2), x], x] + Dist[e/(2*c), Int[(b + 2*c*x)/(a + b*x + c*x^2), x], x] /; FreeQ[{a, b, c, d, e}, x] &
& NeQ[2*c*d - b*e, 0] && NeQ[b^2 - 4*a*c, 0] && !NiceSqrtQ[b^2 - 4*a*c]
Rule 1183
Int[((d_) + (e_.)*(x_)^2)/((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[a/c, 2]}, With[{r =
Rt[2*q - b/c, 2]}, Dist[1/(2*c*q*r), Int[(d*r - (d - e*q)*x)/(q - r*x + x^2), x], x] + Dist[1/(2*c*q*r), Int[(
d*r + (d - e*q)*x)/(q + r*x + x^2), x], x]]] /; FreeQ[{a, b, c, d, e}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2
- b*d*e + a*e^2, 0] && NegQ[b^2 - 4*a*c]
Rule 1435
Int[((d_) + (e_.)*(x_)^(n_))/((a_) + (b_.)*(x_)^(n_) + (c_.)*(x_)^(n2_)), x_Symbol] :> With[{q = Rt[-2*(d/e) -
b/c, 2]}, Dist[e/(2*c*q), Int[(q - 2*x^(n/2))/Simp[d/e + q*x^(n/2) - x^n, x], x], x] + Dist[e/(2*c*q), Int[(q
+ 2*x^(n/2))/Simp[d/e - q*x^(n/2) - x^n, x], x], x]] /; FreeQ[{a, b, c, d, e}, x] && EqQ[n2, 2*n] && NeQ[b^2
- 4*a*c, 0] && EqQ[c*d^2 - a*e^2, 0] && IGtQ[n/2, 0] && !GtQ[b^2 - 4*a*c, 0]
Rubi steps \begin{align*}
\text {integral}& = -\frac {\int \frac {\sqrt {2-b}+2 x^2}{-1-\sqrt {2-b} x^2-x^4} \, dx}{2 \sqrt {2-b}}-\frac {\int \frac {\sqrt {2-b}-2 x^2}{-1+\sqrt {2-b} x^2-x^4} \, dx}{2 \sqrt {2-b}} \\ & = \frac {\int \frac {\sqrt {2-\sqrt {2-b}} \sqrt {2-b}-\left (-2+\sqrt {2-b}\right ) x}{1-\sqrt {2-\sqrt {2-b}} x+x^2} \, dx}{4 \sqrt {2-\sqrt {2-b}} \sqrt {2-b}}+\frac {\int \frac {\sqrt {2-\sqrt {2-b}} \sqrt {2-b}+\left (-2+\sqrt {2-b}\right ) x}{1+\sqrt {2-\sqrt {2-b}} x+x^2} \, dx}{4 \sqrt {2-\sqrt {2-b}} \sqrt {2-b}}+\frac {\int \frac {\sqrt {2+\sqrt {2-b}} \sqrt {2-b}-\left (2+\sqrt {2-b}\right ) x}{1-\sqrt {2+\sqrt {2-b}} x+x^2} \, dx}{4 \sqrt {2+\sqrt {2-b}} \sqrt {2-b}}+\frac {\int \frac {\sqrt {2+\sqrt {2-b}} \sqrt {2-b}+\left (2+\sqrt {2-b}\right ) x}{1+\sqrt {2+\sqrt {2-b}} x+x^2} \, dx}{4 \sqrt {2+\sqrt {2-b}} \sqrt {2-b}} \\ & = -\left (\frac {1}{8} \left (-1+\frac {2}{\sqrt {2-b}}\right ) \int \frac {1}{1-\sqrt {2+\sqrt {2-b}} x+x^2} \, dx\right )-\frac {1}{8} \left (-1+\frac {2}{\sqrt {2-b}}\right ) \int \frac {1}{1+\sqrt {2+\sqrt {2-b}} x+x^2} \, dx+\frac {1}{8} \left (1+\frac {2}{\sqrt {2-b}}\right ) \int \frac {1}{1-\sqrt {2-\sqrt {2-b}} x+x^2} \, dx+\frac {1}{8} \left (1+\frac {2}{\sqrt {2-b}}\right ) \int \frac {1}{1+\sqrt {2-\sqrt {2-b}} x+x^2} \, dx+\frac {\sqrt {2-\sqrt {2-b}} \int \frac {-\sqrt {2-\sqrt {2-b}}+2 x}{1-\sqrt {2-\sqrt {2-b}} x+x^2} \, dx}{8 \sqrt {2-b}}-\frac {\sqrt {2-\sqrt {2-b}} \int \frac {\sqrt {2-\sqrt {2-b}}+2 x}{1+\sqrt {2-\sqrt {2-b}} x+x^2} \, dx}{8 \sqrt {2-b}}-\frac {\sqrt {2+\sqrt {2-b}} \int \frac {-\sqrt {2+\sqrt {2-b}}+2 x}{1-\sqrt {2+\sqrt {2-b}} x+x^2} \, dx}{8 \sqrt {2-b}}+\frac {\sqrt {2+\sqrt {2-b}} \int \frac {\sqrt {2+\sqrt {2-b}}+2 x}{1+\sqrt {2+\sqrt {2-b}} x+x^2} \, dx}{8 \sqrt {2-b}} \\ & = \frac {\sqrt {2-\sqrt {2-b}} \log \left (1-\sqrt {2-\sqrt {2-b}} x+x^2\right )}{8 \sqrt {2-b}}-\frac {\sqrt {2-\sqrt {2-b}} \log \left (1+\sqrt {2-\sqrt {2-b}} x+x^2\right )}{8 \sqrt {2-b}}-\frac {\sqrt {2+\sqrt {2-b}} \log \left (1-\sqrt {2+\sqrt {2-b}} x+x^2\right )}{8 \sqrt {2-b}}+\frac {\sqrt {2+\sqrt {2-b}} \log \left (1+\sqrt {2+\sqrt {2-b}} x+x^2\right )}{8 \sqrt {2-b}}+\frac {1}{4} \left (-1-\frac {2}{\sqrt {2-b}}\right ) \text {Subst}\left (\int \frac {1}{-2-\sqrt {2-b}-x^2} \, dx,x,-\sqrt {2-\sqrt {2-b}}+2 x\right )+\frac {1}{4} \left (-1-\frac {2}{\sqrt {2-b}}\right ) \text {Subst}\left (\int \frac {1}{-2-\sqrt {2-b}-x^2} \, dx,x,\sqrt {2-\sqrt {2-b}}+2 x\right )-\frac {1}{4} \left (1-\frac {2}{\sqrt {2-b}}\right ) \text {Subst}\left (\int \frac {1}{-2+\sqrt {2-b}-x^2} \, dx,x,-\sqrt {2+\sqrt {2-b}}+2 x\right )-\frac {1}{4} \left (1-\frac {2}{\sqrt {2-b}}\right ) \text {Subst}\left (\int \frac {1}{-2+\sqrt {2-b}-x^2} \, dx,x,\sqrt {2+\sqrt {2-b}}+2 x\right ) \\ & = -\frac {\sqrt {2+\sqrt {2-b}} \tan ^{-1}\left (\frac {\sqrt {2-\sqrt {2-b}}-2 x}{\sqrt {2+\sqrt {2-b}}}\right )}{4 \sqrt {2-b}}+\frac {\sqrt {2-\sqrt {2-b}} \tan ^{-1}\left (\frac {\sqrt {2+\sqrt {2-b}}-2 x}{\sqrt {2-\sqrt {2-b}}}\right )}{4 \sqrt {2-b}}+\frac {\sqrt {2+\sqrt {2-b}} \tan ^{-1}\left (\frac {\sqrt {2-\sqrt {2-b}}+2 x}{\sqrt {2+\sqrt {2-b}}}\right )}{4 \sqrt {2-b}}-\frac {\sqrt {2-\sqrt {2-b}} \tan ^{-1}\left (\frac {\sqrt {2+\sqrt {2-b}}+2 x}{\sqrt {2-\sqrt {2-b}}}\right )}{4 \sqrt {2-b}}+\frac {\sqrt {2-\sqrt {2-b}} \log \left (1-\sqrt {2-\sqrt {2-b}} x+x^2\right )}{8 \sqrt {2-b}}-\frac {\sqrt {2-\sqrt {2-b}} \log \left (1+\sqrt {2-\sqrt {2-b}} x+x^2\right )}{8 \sqrt {2-b}}-\frac {\sqrt {2+\sqrt {2-b}} \log \left (1-\sqrt {2+\sqrt {2-b}} x+x^2\right )}{8 \sqrt {2-b}}+\frac {\sqrt {2+\sqrt {2-b}} \log \left (1+\sqrt {2+\sqrt {2-b}} x+x^2\right )}{8 \sqrt {2-b}} \\
\end{align*}
Mathematica [C] (verified)
Result contains higher order function than in optimal. Order 9 vs. order 3 in optimal.
Time = 0.03 (sec) , antiderivative size = 57, normalized size of antiderivative = 0.11
\[
\int \frac {1-x^4}{1+b x^4+x^8} \, dx=-\frac {1}{4} \text {RootSum}\left [1+b \text {$\#$1}^4+\text {$\#$1}^8\&,\frac {-\log (x-\text {$\#$1})+\log (x-\text {$\#$1}) \text {$\#$1}^4}{b \text {$\#$1}^3+2 \text {$\#$1}^7}\&\right ]
\]
[In]
Integrate[(1 - x^4)/(1 + b*x^4 + x^8),x]
[Out]
-1/4*RootSum[1 + b*#1^4 + #1^8 & , (-Log[x - #1] + Log[x - #1]*#1^4)/(b*#1^3 + 2*#1^7) & ]
Maple [C] (verified)
Result contains higher order function than in optimal. Order 9 vs. order 3.
Time = 0.05 (sec) , antiderivative size = 44, normalized size of antiderivative = 0.09
| | |
method | result | size |
| | |
default |
\(\frac {\left (\munderset {\textit {\_R} =\operatorname {RootOf}\left (\textit {\_Z}^{8}+\textit {\_Z}^{4} b +1\right )}{\sum }\frac {\left (-\textit {\_R}^{4}+1\right ) \ln \left (x -\textit {\_R} \right )}{2 \textit {\_R}^{7}+\textit {\_R}^{3} b}\right )}{4}\) |
\(44\) |
risch |
\(\frac {\left (\munderset {\textit {\_R} =\operatorname {RootOf}\left (\textit {\_Z}^{8}+\textit {\_Z}^{4} b +1\right )}{\sum }\frac {\left (-\textit {\_R}^{4}+1\right ) \ln \left (x -\textit {\_R} \right )}{2 \textit {\_R}^{7}+\textit {\_R}^{3} b}\right )}{4}\) |
\(44\) |
| | |
|
|
|
[In]
int((-x^4+1)/(x^8+b*x^4+1),x,method=_RETURNVERBOSE)
[Out]
1/4*sum((-_R^4+1)/(2*_R^7+_R^3*b)*ln(x-_R),_R=RootOf(_Z^8+_Z^4*b+1))
Fricas [B] (verification not implemented)
Leaf count of result is larger than twice the leaf count of optimal. 1177 vs. \(2 (397) = 794\).
Time = 0.28 (sec) , antiderivative size = 1177, normalized size of antiderivative = 2.30
\[
\int \frac {1-x^4}{1+b x^4+x^8} \, dx=\text {Too large to display}
\]
[In]
integrate((-x^4+1)/(x^8+b*x^4+1),x, algorithm="fricas")
[Out]
1/4*sqrt(sqrt(1/2)*sqrt(-((b^2 - 4*b + 4)*sqrt((b + 2)/(b^3 - 6*b^2 + 12*b - 8)) + b)/(b^2 - 4*b + 4)))*log(1/
2*((b^2 - 4*b + 4)*sqrt((b + 2)/(b^3 - 6*b^2 + 12*b - 8)) - b + 2)*sqrt(sqrt(1/2)*sqrt(-((b^2 - 4*b + 4)*sqrt(
(b + 2)/(b^3 - 6*b^2 + 12*b - 8)) + b)/(b^2 - 4*b + 4))) + x) - 1/4*sqrt(sqrt(1/2)*sqrt(-((b^2 - 4*b + 4)*sqrt
((b + 2)/(b^3 - 6*b^2 + 12*b - 8)) + b)/(b^2 - 4*b + 4)))*log(-1/2*((b^2 - 4*b + 4)*sqrt((b + 2)/(b^3 - 6*b^2
+ 12*b - 8)) - b + 2)*sqrt(sqrt(1/2)*sqrt(-((b^2 - 4*b + 4)*sqrt((b + 2)/(b^3 - 6*b^2 + 12*b - 8)) + b)/(b^2 -
4*b + 4))) + x) + 1/4*sqrt(-sqrt(1/2)*sqrt(-((b^2 - 4*b + 4)*sqrt((b + 2)/(b^3 - 6*b^2 + 12*b - 8)) + b)/(b^2
- 4*b + 4)))*log(1/2*((b^2 - 4*b + 4)*sqrt((b + 2)/(b^3 - 6*b^2 + 12*b - 8)) - b + 2)*sqrt(-sqrt(1/2)*sqrt(-(
(b^2 - 4*b + 4)*sqrt((b + 2)/(b^3 - 6*b^2 + 12*b - 8)) + b)/(b^2 - 4*b + 4))) + x) - 1/4*sqrt(-sqrt(1/2)*sqrt(
-((b^2 - 4*b + 4)*sqrt((b + 2)/(b^3 - 6*b^2 + 12*b - 8)) + b)/(b^2 - 4*b + 4)))*log(-1/2*((b^2 - 4*b + 4)*sqrt
((b + 2)/(b^3 - 6*b^2 + 12*b - 8)) - b + 2)*sqrt(-sqrt(1/2)*sqrt(-((b^2 - 4*b + 4)*sqrt((b + 2)/(b^3 - 6*b^2 +
12*b - 8)) + b)/(b^2 - 4*b + 4))) + x) - 1/4*sqrt(sqrt(1/2)*sqrt(((b^2 - 4*b + 4)*sqrt((b + 2)/(b^3 - 6*b^2 +
12*b - 8)) - b)/(b^2 - 4*b + 4)))*log(1/2*((b^2 - 4*b + 4)*sqrt((b + 2)/(b^3 - 6*b^2 + 12*b - 8)) + b - 2)*sq
rt(sqrt(1/2)*sqrt(((b^2 - 4*b + 4)*sqrt((b + 2)/(b^3 - 6*b^2 + 12*b - 8)) - b)/(b^2 - 4*b + 4))) + x) + 1/4*sq
rt(sqrt(1/2)*sqrt(((b^2 - 4*b + 4)*sqrt((b + 2)/(b^3 - 6*b^2 + 12*b - 8)) - b)/(b^2 - 4*b + 4)))*log(-1/2*((b^
2 - 4*b + 4)*sqrt((b + 2)/(b^3 - 6*b^2 + 12*b - 8)) + b - 2)*sqrt(sqrt(1/2)*sqrt(((b^2 - 4*b + 4)*sqrt((b + 2)
/(b^3 - 6*b^2 + 12*b - 8)) - b)/(b^2 - 4*b + 4))) + x) - 1/4*sqrt(-sqrt(1/2)*sqrt(((b^2 - 4*b + 4)*sqrt((b + 2
)/(b^3 - 6*b^2 + 12*b - 8)) - b)/(b^2 - 4*b + 4)))*log(1/2*((b^2 - 4*b + 4)*sqrt((b + 2)/(b^3 - 6*b^2 + 12*b -
8)) + b - 2)*sqrt(-sqrt(1/2)*sqrt(((b^2 - 4*b + 4)*sqrt((b + 2)/(b^3 - 6*b^2 + 12*b - 8)) - b)/(b^2 - 4*b + 4
))) + x) + 1/4*sqrt(-sqrt(1/2)*sqrt(((b^2 - 4*b + 4)*sqrt((b + 2)/(b^3 - 6*b^2 + 12*b - 8)) - b)/(b^2 - 4*b +
4)))*log(-1/2*((b^2 - 4*b + 4)*sqrt((b + 2)/(b^3 - 6*b^2 + 12*b - 8)) + b - 2)*sqrt(-sqrt(1/2)*sqrt(((b^2 - 4*
b + 4)*sqrt((b + 2)/(b^3 - 6*b^2 + 12*b - 8)) - b)/(b^2 - 4*b + 4))) + x)
Sympy [A] (verification not implemented)
Time = 1.87 (sec) , antiderivative size = 76, normalized size of antiderivative = 0.15
\[
\int \frac {1-x^4}{1+b x^4+x^8} \, dx=- \operatorname {RootSum} {\left (t^{8} \cdot \left (65536 b^{4} - 524288 b^{3} + 1572864 b^{2} - 2097152 b + 1048576\right ) + t^{4} \cdot \left (256 b^{3} - 1024 b^{2} + 1024 b\right ) + 1, \left ( t \mapsto t \log {\left (1024 t^{5} b^{2} - 4096 t^{5} b + 4096 t^{5} + 4 t b - 4 t + x \right )} \right )\right )}
\]
[In]
integrate((-x**4+1)/(x**8+b*x**4+1),x)
[Out]
-RootSum(_t**8*(65536*b**4 - 524288*b**3 + 1572864*b**2 - 2097152*b + 1048576) + _t**4*(256*b**3 - 1024*b**2 +
1024*b) + 1, Lambda(_t, _t*log(1024*_t**5*b**2 - 4096*_t**5*b + 4096*_t**5 + 4*_t*b - 4*_t + x)))
Maxima [F]
\[
\int \frac {1-x^4}{1+b x^4+x^8} \, dx=\int { -\frac {x^{4} - 1}{x^{8} + b x^{4} + 1} \,d x }
\]
[In]
integrate((-x^4+1)/(x^8+b*x^4+1),x, algorithm="maxima")
[Out]
-integrate((x^4 - 1)/(x^8 + b*x^4 + 1), x)
Giac [F]
\[
\int \frac {1-x^4}{1+b x^4+x^8} \, dx=\int { -\frac {x^{4} - 1}{x^{8} + b x^{4} + 1} \,d x }
\]
[In]
integrate((-x^4+1)/(x^8+b*x^4+1),x, algorithm="giac")
[Out]
integrate(-(x^4 - 1)/(x^8 + b*x^4 + 1), x)
Mupad [B] (verification not implemented)
Time = 9.53 (sec) , antiderivative size = 5341, normalized size of antiderivative = 10.45
\[
\int \frac {1-x^4}{1+b x^4+x^8} \, dx=\text {Too large to display}
\]
[In]
int(-(x^4 - 1)/(b*x^4 + x^8 + 1),x)
[Out]
- atan((((-(4*b + ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(25
6*b + ((-(4*b + ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(2621
44*b - 196608*b^2 - 196608*b^3 + 49152*b^4 + 49152*b^5 - 4096*b^6 - 4096*b^7 + 262144) + x*(32768*b - 65536*b^
2 - 32768*b^3 + 20480*b^4 + 10240*b^5 - 2048*b^6 - 1024*b^7 + 65536))*(-(4*b + ((b - 2)^5*(b + 2))^(1/2) - 4*b
^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(3/4) - 64*b^3 - 16*b^4 + 256) - x*(32*b + 48*b^2 + 24*b^3
+ 4*b^4))*(-(4*b + ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*1
i - ((-(4*b + ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(256*b
+ ((-(4*b + ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(262144*b
- 196608*b^2 - 196608*b^3 + 49152*b^4 + 49152*b^5 - 4096*b^6 - 4096*b^7 + 262144) - x*(32768*b - 65536*b^2 -
32768*b^3 + 20480*b^4 + 10240*b^5 - 2048*b^6 - 1024*b^7 + 65536))*(-(4*b + ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 +
b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(3/4) - 64*b^3 - 16*b^4 + 256) + x*(32*b + 48*b^2 + 24*b^3 + 4
*b^4))*(-(4*b + ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*1i)/(
((-(4*b + ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(256*b + ((
-(4*b + ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(262144*b - 1
96608*b^2 - 196608*b^3 + 49152*b^4 + 49152*b^5 - 4096*b^6 - 4096*b^7 + 262144) + x*(32768*b - 65536*b^2 - 3276
8*b^3 + 20480*b^4 + 10240*b^5 - 2048*b^6 - 1024*b^7 + 65536))*(-(4*b + ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3
)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(3/4) - 64*b^3 - 16*b^4 + 256) - x*(32*b + 48*b^2 + 24*b^3 + 4*b^4
))*(-(4*b + ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4) + ((-(4*b
+ ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(256*b + ((-(4*b +
((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(262144*b - 196608*b
^2 - 196608*b^3 + 49152*b^4 + 49152*b^5 - 4096*b^6 - 4096*b^7 + 262144) - x*(32768*b - 65536*b^2 - 32768*b^3 +
20480*b^4 + 10240*b^5 - 2048*b^6 - 1024*b^7 + 65536))*(-(4*b + ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*
(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(3/4) - 64*b^3 - 16*b^4 + 256) + x*(32*b + 48*b^2 + 24*b^3 + 4*b^4))*(-(4
*b + ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)))*(-(4*b + ((b -
2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*2i - 2*atan((((-(4*b + ((b
- 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(((-(4*b + ((b - 2)^5*(b
+ 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(262144*b - 196608*b^2 - 196608*b^
3 + 49152*b^4 + 49152*b^5 - 4096*b^6 - 4096*b^7 + 262144)*1i + x*(32768*b - 65536*b^2 - 32768*b^3 + 20480*b^4
+ 10240*b^5 - 2048*b^6 - 1024*b^7 + 65536))*(-(4*b + ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 3
2*b - 8*b^3 + b^4 + 16)))^(3/4)*1i - 256*b + 64*b^3 + 16*b^4 - 256)*1i + x*(32*b + 48*b^2 + 24*b^3 + 4*b^4))*(
-(4*b + ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4) - ((-(4*b + (
(b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(((-(4*b + ((b - 2)^5*
(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(262144*b - 196608*b^2 - 196608*
b^3 + 49152*b^4 + 49152*b^5 - 4096*b^6 - 4096*b^7 + 262144)*1i - x*(32768*b - 65536*b^2 - 32768*b^3 + 20480*b^
4 + 10240*b^5 - 2048*b^6 - 1024*b^7 + 65536))*(-(4*b + ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 -
32*b - 8*b^3 + b^4 + 16)))^(3/4)*1i - 256*b + 64*b^3 + 16*b^4 - 256)*1i - x*(32*b + 48*b^2 + 24*b^3 + 4*b^4))
*(-(4*b + ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4))/(((-(4*b +
((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(((-(4*b + ((b - 2)^
5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(262144*b - 196608*b^2 - 19660
8*b^3 + 49152*b^4 + 49152*b^5 - 4096*b^6 - 4096*b^7 + 262144)*1i + x*(32768*b - 65536*b^2 - 32768*b^3 + 20480*
b^4 + 10240*b^5 - 2048*b^6 - 1024*b^7 + 65536))*(-(4*b + ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2
- 32*b - 8*b^3 + b^4 + 16)))^(3/4)*1i - 256*b + 64*b^3 + 16*b^4 - 256)*1i + x*(32*b + 48*b^2 + 24*b^3 + 4*b^4
))*(-(4*b + ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*1i + ((-(
4*b + ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(((-(4*b + ((b
- 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(262144*b - 196608*b^2 -
196608*b^3 + 49152*b^4 + 49152*b^5 - 4096*b^6 - 4096*b^7 + 262144)*1i - x*(32768*b - 65536*b^2 - 32768*b^3 + 2
0480*b^4 + 10240*b^5 - 2048*b^6 - 1024*b^7 + 65536))*(-(4*b + ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(2
4*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(3/4)*1i - 256*b + 64*b^3 + 16*b^4 - 256)*1i - x*(32*b + 48*b^2 + 24*b^3 +
4*b^4))*(-(4*b + ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*1i))
*(-(4*b + ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4) - atan((((-
(4*b - ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(256*b + ((-(4
*b - ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(262144*b - 1966
08*b^2 - 196608*b^3 + 49152*b^4 + 49152*b^5 - 4096*b^6 - 4096*b^7 + 262144) + x*(32768*b - 65536*b^2 - 32768*b
^3 + 20480*b^4 + 10240*b^5 - 2048*b^6 - 1024*b^7 + 65536))*(-(4*b - ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(
512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(3/4) - 64*b^3 - 16*b^4 + 256) - x*(32*b + 48*b^2 + 24*b^3 + 4*b^4))*
(-(4*b - ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*1i - ((-(4*b
- ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(256*b + ((-(4*b -
((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(262144*b - 196608*b
^2 - 196608*b^3 + 49152*b^4 + 49152*b^5 - 4096*b^6 - 4096*b^7 + 262144) - x*(32768*b - 65536*b^2 - 32768*b^3 +
20480*b^4 + 10240*b^5 - 2048*b^6 - 1024*b^7 + 65536))*(-(4*b - ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*
(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(3/4) - 64*b^3 - 16*b^4 + 256) + x*(32*b + 48*b^2 + 24*b^3 + 4*b^4))*(-(4
*b - ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*1i)/(((-(4*b - (
(b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(256*b + ((-(4*b - ((b
- 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(262144*b - 196608*b^2 -
196608*b^3 + 49152*b^4 + 49152*b^5 - 4096*b^6 - 4096*b^7 + 262144) + x*(32768*b - 65536*b^2 - 32768*b^3 + 204
80*b^4 + 10240*b^5 - 2048*b^6 - 1024*b^7 + 65536))*(-(4*b - ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*
b^2 - 32*b - 8*b^3 + b^4 + 16)))^(3/4) - 64*b^3 - 16*b^4 + 256) - x*(32*b + 48*b^2 + 24*b^3 + 4*b^4))*(-(4*b -
((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4) + ((-(4*b - ((b - 2)
^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(256*b + ((-(4*b - ((b - 2)^5
*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(262144*b - 196608*b^2 - 196608
*b^3 + 49152*b^4 + 49152*b^5 - 4096*b^6 - 4096*b^7 + 262144) - x*(32768*b - 65536*b^2 - 32768*b^3 + 20480*b^4
+ 10240*b^5 - 2048*b^6 - 1024*b^7 + 65536))*(-(4*b - ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 3
2*b - 8*b^3 + b^4 + 16)))^(3/4) - 64*b^3 - 16*b^4 + 256) + x*(32*b + 48*b^2 + 24*b^3 + 4*b^4))*(-(4*b - ((b -
2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)))*(-(4*b - ((b - 2)^5*(b +
2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*2i - 2*atan((((-(4*b - ((b - 2)^5*(b
+ 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(((-(4*b - ((b - 2)^5*(b + 2))^(1/2
) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(262144*b - 196608*b^2 - 196608*b^3 + 49152*b
^4 + 49152*b^5 - 4096*b^6 - 4096*b^7 + 262144)*1i + x*(32768*b - 65536*b^2 - 32768*b^3 + 20480*b^4 + 10240*b^5
- 2048*b^6 - 1024*b^7 + 65536))*(-(4*b - ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3
+ b^4 + 16)))^(3/4)*1i - 256*b + 64*b^3 + 16*b^4 - 256)*1i + x*(32*b + 48*b^2 + 24*b^3 + 4*b^4))*(-(4*b - ((b
- 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4) - ((-(4*b - ((b - 2)^5*(
b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(((-(4*b - ((b - 2)^5*(b + 2))^(1
/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(262144*b - 196608*b^2 - 196608*b^3 + 49152
*b^4 + 49152*b^5 - 4096*b^6 - 4096*b^7 + 262144)*1i - x*(32768*b - 65536*b^2 - 32768*b^3 + 20480*b^4 + 10240*b
^5 - 2048*b^6 - 1024*b^7 + 65536))*(-(4*b - ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b
^3 + b^4 + 16)))^(3/4)*1i - 256*b + 64*b^3 + 16*b^4 - 256)*1i - x*(32*b + 48*b^2 + 24*b^3 + 4*b^4))*(-(4*b - (
(b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4))/(((-(4*b - ((b - 2)^5
*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(((-(4*b - ((b - 2)^5*(b + 2))^
(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(262144*b - 196608*b^2 - 196608*b^3 + 491
52*b^4 + 49152*b^5 - 4096*b^6 - 4096*b^7 + 262144)*1i + x*(32768*b - 65536*b^2 - 32768*b^3 + 20480*b^4 + 10240
*b^5 - 2048*b^6 - 1024*b^7 + 65536))*(-(4*b - ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8
*b^3 + b^4 + 16)))^(3/4)*1i - 256*b + 64*b^3 + 16*b^4 - 256)*1i + x*(32*b + 48*b^2 + 24*b^3 + 4*b^4))*(-(4*b -
((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*1i + ((-(4*b - ((b -
2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(((-(4*b - ((b - 2)^5*(b +
2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*(262144*b - 196608*b^2 - 196608*b^3
+ 49152*b^4 + 49152*b^5 - 4096*b^6 - 4096*b^7 + 262144)*1i - x*(32768*b - 65536*b^2 - 32768*b^3 + 20480*b^4 +
10240*b^5 - 2048*b^6 - 1024*b^7 + 65536))*(-(4*b - ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*
b - 8*b^3 + b^4 + 16)))^(3/4)*1i - 256*b + 64*b^3 + 16*b^4 - 256)*1i - x*(32*b + 48*b^2 + 24*b^3 + 4*b^4))*(-(
4*b - ((b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)*1i))*(-(4*b - (
(b - 2)^5*(b + 2))^(1/2) - 4*b^2 + b^3)/(512*(24*b^2 - 32*b - 8*b^3 + b^4 + 16)))^(1/4)